One-dimensional Josephson arrays as superlattices for single Cooper pairs 
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We investigate uniform one-dimensional arrays of small Josephson junctions (Ej <C Ec, 
Ec = (2e) 2 /2C) with a realistic Coulomb interaction U(x) — EcXexp(— \x\/X) (here A 2> 1 is 
the screening length in units of the lattice constant of the array). At low energies this system 
can be described in terms of interacting Bose particles (extra single Cooper pairs) on the lattice. 
With increasing concentration v of extra Cooper pairs, a crossover from the Bose gas phase to 
the Wigner crystal phase and then to the superlattice regime occurs. The phase diagram in the 
superlattice regime consists of commensurable insulating phases with v = 1/1 (I is integer) sepa- 
rated by superconducting regions where the current is carried by excitations with fractional electric 
charge q — ±2e/i. The Josephson current through a ring-shaped array pierced by magnetic flux is 
calculated for all of the phases. 

PACS numbers: 05.30.Jp, 73.40.Gk, 74.50.+r, 73.20.Dx, 72.15.Nj 



The investigation of arrays of small Josephson junc- 
tions attracts growing interest of theoreticians and ex- 
perimentalists (see Q for a review). In such arrays the 
Coulomb energy Ec — (2e) 2 /2C can be of the order of 
the Josephson energy Ej (here C is the capacitance of 
the junctions). Since the Josephson phase <f> and the elec- 
tric charge Q on the islands are canonically conjugated 
operators, [<f>, Q] = 2ei, a number of nontrivial quantum 
phenomena arise due to a competition between the phase 
(or vortex) and charge degrees of freedom. In particular, 
the point of the superconductor-insulator transition de- 
pends on the magnetic frustration B a nd on the electro- 
chemical potential fj, of the array fflq]. 

Existing theories of electric field-induced supercon- 
ductor-insulator transition |3]-[3| treat predominantly the 
cases of on-site or nearest-neighbor Coulomb interaction, 
although for typical experimental parameters the range 
of the interaction is large JjJ , A = 3 -j- 20 (in units of the 
lattice constant of the array) . The phase diagram of the 
superconductor-insulator transition becomes rather com- 
plicated for A 3> 1, and large concentrations v ~ 1/A d of 
field-induced ("extra") Cooper pairs H (here d is a di- 
mension of the array). This is related to the fact that the 
extra Cooper pairs (ECP) can form a variety of configu- 
rations commensurable with the lattice of junctions. On 
the other hand, it is known that at very small con- 
centrations |P v <C 1/A d the ECP form a Bose gas with 
hard-core interaction and commensurability with the lat- 
tice plays no role. In this work we study the transition 
between these two very different regimes. 

We consider uniform one-dimensional (1-D) Josephson 
arrays with large Coulomb energy Ec 3> Ej and long 
range Coulomb interactions, A » 1. We focus on the 
regime with low concentrations v <C 1/A of ECP. We 
will show that the phase diagram of the superconductor- 
insulator transition has simple structure in this regime. 
In the insulating phase the ECP form a regular superlat- 



tice with the period I (v — I /I). The system starts show- 
ing superconducting properties when the first mobile ex- 
citation appears. The latter corresponds to a change of 
the distance between two neighboring ECP in a super- 
lattice by Al — ±1. This excitation can be treated as a 
quasiparticle with fractional charge q — ^f2e/l. The su- 
perconducting phase can be viewed a gas of such quasi- 
particles on the lattice. 

The Hamiltonian H = He + Hj of 1-D Josephson 
array consists of a Coulomb term He and a Joseph- 
son term Hj. The Coulomb energy is given by He = 
\ Si j=i niU(i—j)nj — fiN, where ni is a (positive or neg- 
ative) number of ECP on the electrode i, N — Yli=i n i> 
and U(x) = Ec Xcxp(— x/X) is the Coulomb interaction 
between ECP. The interaction is screened at the length 
A = (C/Cg) 1 / 2 due to a finite self-capacitance C g of the 
islands. The gate voltage V plays the role of the chemical 
potential, fi — 2eV. 

We consider first the limit of zero Josephson coupling. 
For \fi\ < ix tr = XEc/2 there is no ECP in the array 
(ni = for all i). Just above the threshold, < [i~[itr *C 
fitr, the ground state is still characterized by the absence 
of ECP on most of the islands (ni = 0). The rest of 
the islands are occupied by one ECP (ni = 1). The 
configurations with n$ ^ 0, 1 (for some i) are separated 
from the ground state by a Coulomb gap Ac = Ec/X, 
which corresponds to the difference between the Coulomb 
energies of the configurations {n^ = (0, 1, —1, 1, 0, ...) 
and (0, 1, 0, ...). In what follows we restrict the space 
of states to low energy configurations with ni = or 
1. These configurations can be fully characterized by 
the coordinates Xj of ECP on the lattice (xj are integer 
numbers). 

The Josephson term of the Hamiltonian has a stan- 
dard form, Hj = —Ej^2f =1 cos((f>i+i — <f>i — a). Here 
cf>j is the operator of the Josephson phase of the island 
j obeying the commutation relation [nj,4>ji] ~ —Sj,j>i, 
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and a = (27r/L)$/<&o is the vector potential (for circular 
array pierced by a magnetic flux <E>). The Josephson term 
acting on the restricted space of states describes a hop- 
ing of ECP on the lattice with the amplitude Eje ±la /2. 
Corrections to the tunneling amplitude due to the states 
with n ^ 0,1 are small for Ej <C Ac. Therefore, the 
original Hamiltonian H in the low-energy space is equiv- 
alent to the Hamiltonian of Bose particles (ECP) on the 
lattice M, 



N 



H = -Ej cos(pj 
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Xj )-pN, (1) 



where pj are quasimomenta of ECP and p = /j, — fx tr . 
Since the interaction potential U (x) is screened strongly 
on the interparticle distances v~ x 3> A we take into ac- 
count the interaction of the neighboring ECP only. 

Let's consider now the ground state of the system as 
a function of the chemical potential p at fixed Ej. The 
first ECP enters the array at p = —Ej We denote the 
deviation from this threshold by \x = p + Ej. At low 
concentrations v (see Eq.(|])) ECP form a Bose gas with 
hard core interaction H (Fig. [I]). The chemical potential 
p,i BG ) — (tt 2 /2)Ejv 2 is determined by the kinetic energy 
of the particles in parabolic band [fL0|. Taking into ac- 
count a finite range of the interaction we obtain by a 
variational calculation that the correction to 

p(BG) due 

to a finite size d = Aln(A 3 Ec / Ej) of the core is small for 



v\]n(\ 3 E c /Ej) < 1. 



(2) 



With increasing concentration ECP get localized 
in coordinate space. We assume apriori that 
ECP form 1-D Wigner crystal. Expanding ki- 
netic and potential energy (^) up to quadratic 
terms we obtain the chemical potential p^ wc ^ = 
(2/ir)(E c Ej/X) 1 / 2 exp(-l/2Xv) + XE C exp{-l/ Xv) and 
the mean-square displacement of neighboring particles 
((Xi+i - Xl ) 2 ) = ^Ej/U^- 1 )] 1 / 2 /^, where 

U (l) = (l/2)d 2 U(x)/dx 2 \ x= i = (£ c /2A)exp(-//A). (3) 

Note that the concentration of ECP v(p) increases much 
more slowly than in the Bose gas regime. This expan- 
sion is legitimate if the fluctuations of displacement are 
small, {{x i+ i — Xi) 2 ) <C A 2 , and the kinetic energy per 
ECP is much less than Ej. These conditions determine 
respectively the lower and the upper bound of the range 
of concentrations 
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in which the Wigner crystal phase exists (Fig. |l]). For 
both the phases (Bose gas and Wigner crystal) the Hamil- 
tonian is quadratic in momenta of the particles. There- 
fore, the vector potential a is coupled to the momenta 
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State: Bose gas Wigner crystal Superlattice 

Interaction: unharmonic harmonic harmonic 

Band: parabolic parabolic cos p 



FIG. 1. Schematic phase diagram of 1-D Josephson array. 

of the center of mass only. For this reason the Joseph- 
son current through a ring-shaped array is given by a 
universal expression 
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vj = nvEj, 
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(for |$| < $0/2)- A similar result holds for the persistent 
current of interacting Fermions JHJ vj being the Fermi 
velocity in that case. 

For larger concentrations of ECP the kinetic energy 
per particle becomes comparable to the bandwidth Ej. 
To investigate this case we start from the limit Ej — > 0. 
First we determine the range of the chemical potential 
/ij ; _(0) < p < /iz,+ (0) where the configuration with 
equidistantly spaced (at a distance I) ECP is the ground 
state of the system. In order to add (subtract) a Cooper 
pair into (from) this configuration, one has to convert I 
bonds of length I between neighboring ECP into I bonds 
of length I — 1 (I + 1). The energy required for this 
conversion is given by /i;,±(0) = ±l[U(l =F 1) — U(l)] ~ 
e/(l ± 1/2A) with e t = E c le~ l / X . Clearly, addition of 
the second, third, etc. particle to the system requires 
the same energy ^+(0). For this reason, the ground 
state corresponds to a regular superlattice of ECP in the 
array; the period of the superlattice changes abruptly 
from I to / — 1 at p = ^i,+(0) = ^i/_i i _(0), see Fig. ||a. 
This simple picture of the ground state is valid if one 
takes into account the interaction of neighboring ECP 
only (see Eq.(Q)). The interaction of next-nearest neigh- 
bors will lead to new ground states |12j| in exponentially 
narrow regions \p — [ii t +(0)\ ~ e/e _ '' A near the points 
R+(0). 

For small but finite Josephson coupling, Ej <C Uq{1), 
each bond of the length I ± 1 can be considered as a mo- 
bile excitation above the ground state where all the bonds 
have the length I. Since a shift of an ECP to the neigh- 
boring island of the array leads to a shift of the excitation 
by =fl lattice cells, the excitations have fractional charge 
=f2e/l ]l3| . Tunneling of excitations decreases the energy 
by an amount Ej per excitation. Since addition (subtrac- 
tion) of one Cooper pair to (from) the array is accompa- 
nied by the creation of I excitations, the threshold chem- 
ical potentials are given by ni,±(Ej) — ///,±(0) =F lEj. 

In the superlattice regime under consideration (see 
Fig. |l|) the fluctuations of displacement of ECP are 
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small and we can again expand the interaction term of 
the Hamiltonian ([[]) around the average distance I be- 
tween the particles. We introduce the (integer) deviation 
— I of the distance between neighboring 
ECP from I, and the canonically conjugated operator ifj, 
so that pj = ifj+i — ifj. The number of particles N can 
be expressed via nj, N = (L — ^2 rij)/l. In terms of the 
new variables the Hamiltonian (|l|) can be written as 

H = -Ej^2coa(tpj+i - <Pj - a) + C/oW^K + <^) 2 , 
3 i 

(6) 

where Sfi = A(/2— e;)/e;. This Hamiltonian formally coin- 
cides with the Hamiltonian of 1-D Josephson array with 
on-site Coulomb interactions, which has been extensively 
studied in Refs. §-§@. 

The boundary of the commensurable phase can be de- 
termined from the comparison of the ground state en- 
ergies of the Hamiltonian (^|) in two subspaces of states 
with Y^ n j = ano - with = ±1 (the sign coincides 

with the sign of e; — p, in (p)). In the limit Ej — > 
the ground states in these subspaces are given by = 
|0, 0, ...) and * ±1 = 2V- 1 /2£ i |0,...,nj = ±1,0,...). 
Evaluating the energies of these states up to the third 
order in Ej (see Ref. M) we obtain 
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for Ej <C Uq(1), see Fig. ||d. Note that the term linear 
in Ej coincides with that obtained above from a simple 
argument. 

With increasing Ej the range of the chemical po- 
tential corresponding to commensurable phase decreases 
and both phase boundaries tend to the critical point, 
^i,±{Ej Cr ^) — » ei, see Fig. ^b. Clearly, the true behav- 
ior near the critical point cannot be described by per- 
turbation theory of finite order. To extent the pertur- 
bative approach, an extrapolation to infinite order in 
Ej was proposed 0. Unfortunately this (somewhat 
speculative) extrapolation fails to converge to a criti- 
cal point for the 1-D system. To determine the be- 
haviour near the critical point E^ r ^ one can map the 
Hamiltonian @ (with Ji = ej) onto 2-D XY model Q. 
The parameter (2Uo(l)/Ej) 1 ^ 2 plays a role of dimension- 
less temperature ksT/J in the XY model. The point 
of the Kosterlitz-Thouless transition pl| corresponds to 
see Fig. |2|b. Below the transition tem- 
perature (Ej > E^ r ^) spin-spin correlations in the XY 
model decay algebraically with distance. The Josephson 
array shows superconducting properties: the Josephson 
current is inversely proportional to L. It scales as 

E { j cr) + (here c is 
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FIG. 2. Phase diagram of the superconductor-insulator 
transition in 1-D Josephson array (A = 3). (a) The over- 
all view. The insulating spike-like regions from right to left 
correspond to v — 1/7, 1/8, .... (b) Boundary of the insulat- 
ing phase with v = l/l (I 1). The curves from bottom 
to top correspond to the results of the mean-field approach, 
the third order perturbation theory, and the extrapolation 
of the perturbation theory to infinite order connected with 
the extrapolation of the Kosterlitz-Thouless scaling, (c) Con- 
centration v(fi) of ECP and critical Josephson current I c (fb) 
in the superconducting region 1/8 < v < 1/7 (I = 8). We 
normalize I c by the critical current (5) for the concentration 
v = l/l of ECP, I C}1 = 2neEj/lL. 

universal constant). Above the transition temperature 
(Ej < E^ ) the correlations in the XY model decay ex- 



Ey>-0) the 



exp{-b[E J cr) /(E J cr) 



poncntially. Near the critical point ( Ej 
coherence length is given by £ 

Ej)] 1 / 2 }, where b ~ 2. As a result, the Josephson cur- 
rent through 1-D array decays as exp(— L/l£) signalling 
the formation of the insulating phase. 

In the insulating phase, the energy gap for the mo- 
bile excitations scales as J3|,^| For this reason, the 
boundary of insulating phase near the critical point is 
given by 



l H, ± (E J )=e l ±alU (l)/^ 



(8) 



1 



'Ej-E [ T cr) at E, 



with a ~ 1. The upper curve in Fig. ^> corresponds to an 
extrapolation of Eq. (|^) from the neighborhood of the crit- 
non- ical point to lower values of Ej (1.56C/ (0 < Ej < Ej Cr ^). 
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This extrapolation is joined to the extrapolation of the 
perturbative result (0) to infinite order in Ej (drawn in 
the range Ej < 1.56J7o(Z)). A smooth connection of the 
two curves occurs for a ~ 1.73 (and 6 = 2). 

For completeness, we present the mean-field result for 
the phase boundary 



M,±(£j) 




(9) 



Although this expression coincides with Eq. ^ to first or- 
der in Ej, the overall shape of the boundary (Fig. ||b) 
and the critical value of the Josephson energy E^ dif- 
fer considerably from the results discussed above. The 
reason for the failure of the mean-field approach is the 
absence of long-range order in one dimension. 

We return now to the consideration of the super- 
conducting phase and concentrate on the case of small 
Josephson coupling, Ej <C Uq(1). The commensurable 
phases with v = l/l and v = 1/(1 — 1) are separated 
by a narrow superconducting region, /j,i t +(Ej) < [i < 
Hl-\-(Ej), of the width (21 — l)Ej. In this region the 
ground state of the array can be viewed as a gas of frac- 
tionally charged mobile excitations (bonds of length I— 1). 
The excitations interact with each other via a contact po- 
tential U exc (x) — 2Uq(1—\)8 x _o- This is effectively a hard- 
core interaction provided that Ej <C Uq(1). From the 
consideration of the hard-core particles (bonds of length 
I — 1) on a lattice formed by the bonds of length I we 
obtain the chemical potential p and the Josephson cur- 
rent I j as functions of the concentration v = 1/(1 — q) of 
ECP, 

p — /Lt;,+ (0) — Ej{(l — q) cosirq + ir^ 1 sin7rg}, (10) 



li 



4eEj sin irq <E> 
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where the parameter q takes on values in the range < 
q < 1. These dependences are presented in Fig. 

Finally, we discuss effects which are specific for finite 
size circular arrays. If the size L of the array is commen- 
surable with the spacing I of the superlattice of ECP, 
the Josephson current is exponentially small in the insu- 
lating phase. However, if L/l is not an integer, a num- 
ber of residual mobile excitations remain in the array 
in the insulating phase, ^l^-(Ej) < p < fii t +(Ej). In 
the lower part of this range, pi t -(Ej) < p < pi t o(Ej), 
there are m = mod(L,l) residual excitations (bonds of 
the length I + 1) in the ground state. One ECP en- 
ters into the array at p = fii (Ej). As a result, for 
Hi q(Ej) < p < (ii i+ (Ej) a new ground state will contain 
I — m residual excitations (bonds of the length I — 1). 
The threshold chemical potential is given by ni^(Ej) = 
[mfii^(Ej) + (I — m)(ii t+ (Ej)]/l. Since each excitation 
contribute an amount I\cp = (47rei?j/?iL 2 )($/$o) to 



the Josephson current (cf. Eq. (0) with v = 1/L), the 
latter shows a jump at p, — fiifi(Ej). 

There is clearly a need for future investigations, such 
as an analysis of effects of disorder due to the offset 
charges (potential disorder) and due to non-uniformity 
of the Josephson coupling (kinetic disorder). 
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